Bi-stable tunneling current through a molecular quantum dot 
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An exact solution is presented for tunneling through a negative-C/ degenerate molecular quantum 
dot weakly coupled to electrical leads. The tunnel current exhibits hysteresis if the level degeneracy 
of the negative-?/ dot is larger than two. Switching occurs in the voltage range Vi < V < V2 as a 
result of attractive electron correlations in the molecule, which open up a new conducting channel 
when the voltage is above the threshold bias voltage V%. Once this current has been established, the 
extra channel remains open as the voltage is reduced down to the lower threshold voltage Vi . Possible 
realizations of bi-stable molecular quantum dots are fullerenes, especially Ceo, and mixed- valence 
compounds. 

PACS: 85.65.+h, 73.63.-b,73.63.Nm, 71.38.Mx 



Molecular-scale electronics is currently a very active 
area of research jij. The present goals for this field are 
to design and characterize molecules that could be the 
"transmission lines" |2|,[| and active elements in elec- 
tronic circuitry . The dominant mechanism of trans- 
port through active devices will most likely be resonant 
tunneling through electronic molecular states Q (see also 
P,p[ and references therein). A few experimental stud- 
ies provide evidence for various molecular switching 
effects, where the current-voltage (TV) characteristics 
show two branches with high and low current for the same 
voltage. This remarkable phenomenon can result from a 
conformational transformation of certain molecules con- 
taining a "moving part" like a bypyridinium ring, which 
changes its position if the voltage is sufficiently high. The 
transformation necessarily involves a large displacement 
of many atoms so that this ionic switching is rather slow, 
perhaps operating on a millisecond scale. Switching has 
also been observed for simple molecules (organic molecu- 
lar films), and in some cases is strongly dependent on the 
choice of contacts and substrates (?]]. Molecular devices 
that exhibit bi-stability and fast switching could be the 
basis of future oscillators, amplifiers and other impor- 
tant circuit elements. Thus, further progress in molec- 
ular electronics will depend upon finding molecules and 
understanding intrinsic mechanisms for their reversible 
switching from low- to high-current states. 

In this Letter, we study a model quantum dot, which 
exhibits an intrinsic electronic switching of the current 
state due to attractive electron correlations. We show 
that if the degeneracy is larger than two, the tunnel cur- 
rent becomes bistable in some voltage range and the dot 
exhibits a current hysteresis as a function of bias volt- 
age. In the simplest case of a doubly degenerate level, 
the bistability does not occur. We present the exact so- 
lution of the model, allowing for a detailed analysis of 
the current bistability. 

Repulsive electron correlations cause the "Coulomb 



blockade" in the TV characteristics of quantum dots M. 
However, they cannot cause any switching. Here, we 
show that a negative Hubbard U of any origin can provide 
an intrinsic non-retarded current switching of a molecu- 
lar quantum dot. One mechanism, that can produce a 
negative U in molecular systems, is a strong electron- 
phonon (vibronic) interaction. If the tunneling time is 
comparable to or larger than the characteristic phonon 
times, a polaron is formed inside the molecular wire ||. 
There is a wide range of bulk molecular conductors with 
polaronic carriers. Since the formation of polarons in 
polyacetylene (PA) was theoretically discussed Jl(]] , they 
were detected optically in PA ]Tlj , in conjugated poly- 
mers such as polyphenylene, polypyrrole, polythiophcnc, 
polyphenylene sulfid e Jl2[ , Cs-doped biphenyl |L3], n- 
doped bithiophene ]14|, polyphenylenevinylene(PPV)- 
based light emitting diodes H|, and other molecular 
systems. In contrast to bare electrons, polarons attract 
each other at short distances of the order of the inter- 
atomic spacing and form small bipolarons Bipo- 
laron formation can strongly affect the transport prop- 
erties of long molecular wires, as discussed recently [0 . 
When bipolarons are not formed in molecular quantum 
dots because of the short life-time of the carrier inside 
the molecule, the attractive correlations between carri- 
ers still remain. Moreover, attractive short-range corre- 
lations (negative Hubbard U) are feasible even without 
electron-phonon interactions. For example, they might 
be of a pure "chemical" origin, as in the mixed valence 
complexes fl8| . 

Our starting point is the tunneling Hamiltonian, which 
includes a negative Hubbard U in the molecular cigen- 
state £a coupled with the left and right leads by the 



hopping integrals h a 
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Here a a k and c are the annihilation operators in the left 
(a = 1) and right (a — 2) leads, and in the molecule, re- 
spectively, rip = cJjCp, Cafe is the energy dispersion in the 
leads, and U < 0. This negative J7 Hamiltonian is derived 
from the general microscopic Hamiltonian with the bare 
electron-phonon (vibron) interaction and Coulomb repul- 
sion using the canonical transformation Jig]. It was suc- 
cessfully applied to glassy semiconductors |19| , p0[ | , high T c 
superconductors 0,|2l| including doped fullerenes J2^] 
mixed valence compounds Jig ], and more recently to 
superconductor-insulator-superconductor tunneling p3[ . 
The current through the molecular quantum dot is con- 
veniently expressed in terms of the molecular density of 
states (DOS) p M (w) as §] 



duj - f 2 (ui)]'^2r tl (w)p tl (w), 



(2) 



where / 1>2 (w) = {exp[(w + A t eV/2)/T] + 
A is the position of the lowest unoccupied molecular 
level with respect to the chemical potential, r M (w) = 
T lfl (aj)T 2Ll (uj)/[T lfl (uj)) + r 2/i (w)], and r aM (w) = 
27r S Q fe |/iQfe M | 2 ^(^ - Cafe)- Molecular DOS is given by 



Pfj,(uj) = — (1/tt) ImG^w), where G^^uj) is the Fourier 
transform of the r-particle retarded Green's function 
(GF) defined as 



0"), 



-<0(t) Yl ' 



(3) 



i=l 



for 2 < r < oo and G^\t) = -i0{t) ({c M (i), ct }> . Here 
{...,...} is the anticommutator, and 9(t) = 1 for t > 
and zero otherwise. In the following, we apply pertur- 
bation theory with respect to the hopping integrals, ne- 
glecting any contribution to the current other than h 2 , 
but keeping all orders of the negative Hubbard U. Terms 
of higher order in h cannot change the gross I-V features 
for any voltage except the narrow transition regime from 
one lead to another. Then, applying the equations of mo- 
tion for the Heisenberg operators c^t), n^(t) and a a k{t), 
we obtain an infinite set of coupled equations for the 
molecular GFs as 



dGt\t) 



r-l 



dt 



*w e rK>(o) 



[e^ + (r^l)U]G^(t) + UG^ +1 \t), (4) 



where n^(t) = (ct(i)c^(t)^ is the expectation number 
of electrons on the conducting molecular level. For the 
sake of analytical transparency we solve this system for a 
molecule having one d-fold degenerate energy level with 
£ = 0. In this case the set appears to be finite, and it can 



be solved by using the Fourier transformation. Fourier 
transforming Eqs.(4) we find the one-particle GF as 



d-1 



Z r (n) 



r=0 



uj — rU + iS ' 



(5) 



where 6 = +0, n = 71^(0), and 

This is an exact solution with respect to correlations 
which satisfies all sum rules. The electron density n^(t) 
obeys the rate equation, which is obtained by using the 
equations of motion as 



dt 



2^WImA^(t), 



(6) 



where 



r-l 



4&(*)= E KWn^KfcW)- (7) 

These correlation functions should be calculated to first 
order with respect to the hopping integrals h. In this 
order, they satisfy the infinite set of coupled equations 



dt 



r-l 

= hak^n^t) - /(Cafe)] E IT 7l ^( t ) + 

UA^\t) + [U - (r - l)U]A^{t). (8) 



But if we have a finite number of molecular states, the 
set is finite like Eqs.(4). One readily solves this set in 
the steady state, when n. w (t) and (t) become time- 
independent. For a ci-fold degenerate energy level, the 
one-particle correlation function is found as 



4feM = [™-/(£afe)]^E 



Z r (n) 



iak - rU + id 



(9) 



Substituting Eq. (||) into Eq. (|[), we obtain the steady 
state equation for the electron density on the molecule as 



d-i 



^2^2r a (rU)[n-f a (rU)]Z r (n)=0. 



(10) 



a r—0 



Here we assume that T afl (uj) = T a (u>) does not depend 
on /i, otherwise the degeneracy would be removed. To 
simplify the mathematics further, we now assume that 
r a (uj) = r is a constant. Then, from Eq. (||) and Eq. (||) 
the current is found as 



d-i 



3=Y}h{rU)- f 2 (rU)]Z r {n), 



(11) 



r=0 



2 



where j = I/Iq with Iq = edT/2. Let us consider two- 
fold, four-fold, and six-fold degenerate molecular level. 
For d = 2 the kinetic equation is linear in n, and there is 
only one solution, 



E«/«(o) 



2 +E a [/«(«) -/«(tor 



(12) 



The current through a two-fold degenerate molecular dot 
is found as 



3 = 2 



/i(o)[i-/ 2 (t/)]-/ 2 (o)[i-/i([/)] 

2 + E«[/«(0)-/a(E0] 



(13) 



There is no current bistability in this case. Moreover, if 
the temperature is low (T <C A, |?7|) there is practically 
no effect of correlations on the current, j w V#(e|y| — 
2A)/\V\. 

Remarkably, four-fold or higher- degenerate negative 
U dots reveal a switching effect. In this case, the kinetic 
equation is nonlinear, allowing for a few solutions. If 
eV < 2(A — \U\), the only physically allowed solution of 
Eq.(lO) for d = 4 and d = 6 at zero temperature is n = 0. 
If 2(A - \U\) < eV < 2A, T = and \U\ < 2A/d, the 
kinetic equation is reduced to 



2n=l-(l-n) d - 1 . 



(14) 



For d = 4 it has two physical roots, n = and n = 
(3 - 5 1 / 2 )/2 « 0.38 H). In this voltage range /i(0) = 
h(rU) = 0, but h(U) = /!(2C/) = /ifSC/) = 1 at T = 0, 
when \U\ < A/2. Using the sum rule Er=o ^r{ n ) = 1 
and the kinetic equation (10), the current is simplified 
in this voltage range as j — 2n. Hence we obtain two 
stationary states of the molecule with low (zero at T = 0) 
and high current, / « 0.76/o for the same voltage in the 
range 2(A - |Z7|) < eV < 2A. For d = 6, the kinetic 
equation has two physical roots in this voltage range, 
n = and n 0.48, which corresponds to / = and 
/ w 0.96/o, respectively. Above the standard threshold, 
eV > 2A, where fi(rll) = 1 and f 2 (rU) = 0, when 
\U\ < A/3, there is only one solution, n = 0.5 with the 
current I = I a . 

One can better understand the origin of the switch- 
ing phenomenon by taking the limit d » 1. The phys- 
ical roots of Eq. ([li]) are n = and n = 0.5 in this 
limit with the current 1 = and I = Iq, respec- 
tively. This is precisely the solution of the problem in 
the mean- field approximation (MFA), which is a reason- 
able approximation for » 1. Indeed, using MFA one 
replaces the exact two-body interaction in the Hamil- 
tonian for a mean-field potential as ^U^2 n^fi^i w 



,, , n n,. 



/ . Then the MFA DOS 



is given by Pfi(^) = 5[uj — U(d — l)n]. Using the Fermi- 
Dirac Golden rule the rate equation for n becomes 



dn 
~dl 



-2Tn + T^f a [nU{d- 1)]. 



(15) 



For T — there are two stationary solutions of Eq. (|lj) , 
n = and n = 0.5 in the voltage range 2(A — \U\) < 
eV < 2A, and only one solution, n = 0.5 for eV > 2A, 
where U = U(d — l)/2. The MFA current is found as 



j = h{2nU) - h{2nU). 



(16) 



Combining this equation and the rate equation (15) with 
dn/dt = , we obtain the I-V characteristic equation as 



|c/| n m-i T 
— (i _ - 1 - - 



x In 
where 



(1 + R) sinh (eV/2T) 



cosh (eV/2T) 



(17) 



R=\[l-j coth(el//2T)] 



1/2 



sinli (eV/2T) 



1/2 



(18) 



The I /V curves are shown in Fig.l for different temper- 
atures and \U\ = 0.9A. Interestingly, the temperature 
narrows the voltage range of the hysteresis loop, but the 
transition from the low (high)-current branch to the high 
(low)-current branch remains discontinuous. 

Let us examine the stability of each branch in the 
framework of the MFA rate equation jl5|). Introduc- 
ing small fluctuations of the electron density as n(t) = 
n + Snexp^t) and linearizing Eq.(15) with respect to dn 
we find the increment 7, 



= -2r 



\u\r 

2T 



\U\T 
2T 



cosh 



cosh 



( A - 2n\U\ - eV/2 

2n\U\ +eU/2\ 
2T j ' 



(19) 



One can see from this equation that at temperatures T 
\U\ the low-current branch (n = 0) looses its stability at 
the threshold V2 = 2A/e, while the high-current branch 
looses its stability at V\ = 2 (A — \U\)/e. In the voltage 
range Vt < V < V% both branches are stable, 7 w — 2r < 
0. 

Finally, let us analyze the effect of a splitting of the 
degenerate molecular level on the bi-stability. The de- 
generacy could be removed because of Jahn- Teller dis- 
tortions and/or the coupling with the leads. We assume 
that d > 1 levels are evenly distributed in a band of a 
width W. Then the MFA rate equation (15) is modified 
as 

^ = -2Tn, + rJ2Ue. + NU), (20) 

a 

where N = Efi> n n' • ^ or T = in the stationary regime 
(dn,j,/dt = 0) it has two solutions, N = and N = d/2, 



3 
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FIG. 1. The current-voltage hysteresis loop for the 
negative-?/ molecular quantum dot for the molecular level 
with the degeneracies d — 4 and d > 1. The two cases 
are very different in some aspects. Top panel: d = 4, 
|C/|/A = 0.4. In the hysteretic region the current jumps to 
the value I/I Q = 0.76 at the lower threshold (eVi/2A = 0.6) 
and then to the value I/Io = 1 at the higher threshold 
(eV2/2A = 1). Bottom panel: highly degenerate molecular 
level with > 1, |t^|/A = 0.4. The current jumps directly 
to the value I/Io = 1 at the lower threshold. Temperature 
quickly reduces the width of the hysteretic region, much faster 
for d — 4 than for d > 1. In the latter case the hysteresis 
is over at about T/A = 0.1. Inset in the top panel shows a 
schematic of the energy level diagram for d = 4 quantum dot 
under bias voltage V. 



in the voltage range V\ + W < V < V% with the current 
j = and j = 1, respectively. We conclude that the level 
splitting W < \U\ leads to a narrowing of the voltage 
range of the bi-stability similar to the temperature nar- 
rowing shown in Fig.l. A parameter- free estimate of the 
negative Hubbard U in some oxides yields \U\ about a 
few tens of eV |25j . We expect a negative U of the same 
order of magnitude in carbon-based compounds. Among 
potential candidates for the negative U quantum dot are 
a single C^o molecule (d = 6), where the electron- vibronic 
coupling proved to be particularly strong p6|, or other 
fullerenes including short nanotubes (cl > 1) connected 
to metal electrodes. Other likely candidates are mixed- 
valence molecular complexes [ f[8| |. There should be no 
retardation of the switching on the time scale above the 
inverse vibron (phonon) frequency, which is 10 _14 s or less 
in carbon-based compounds. 

In conclusion, we have introduced and solved a model 
for tunneling through a negative-?/ degenerate molecu- 
lar dot weakly coupled to electrical leads. The exact 



many-particle Green's functions of a d— fold degenerate 
molecular level, the density of states and the nonlinear 
rate equation for the electron density on the molecule 
have been derived. We have found the exact solutions 
for the carrier population in the dot and the current for 
the degeneracies of the molecular level d — 2,4,6, and 
ii > 1. The current-voltage characteristics show a hys- 
teretic behavior for d > 2 over a finite voltage range. 
When the voltage increases from zero, the molecule re- 
mains in a low-current state until the threshold Va is 
reached. Remarkably, when the voltage decreases from 
the value above the threshold 2A/e, the molecule re- 
mains in the high-current state down to the voltage 
Vi =(2A - \U\)/e, well below the threshold V 2 . This 
mechanism for electronic molecular switching without re- 
tardation requires many-particle attractive correlations, 
which can arise from strong electron-vibronic coupling 
and/or mixed valence states. Experimental verification 
of such bi-stable systems will require careful collection 
and analysis of both forward and reverse voltage sweeps 
of the tunneling current through candidate molecules. 
The forward voltage sweep by itself will resemble a stan- 
dard Coulomb blockade I-V characteristic with a turn-on 
voltage of V2, whereas the reverse sweep should reveal 
hysteresis. 

This work has been partly supported by DARPA. 
The authors acknowledge useful discussions with P. Ko- 
rnilovich and D. Stewart. 
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